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In  a  series  of  papers  dating  from  1955,  we  have  used  the 
theory  of  invariant  imbedding  to  derive  a  variety  of  nonlinear 
functional  equations  for  the  description  of  radiative  transfer, 
neutron  transport,  diffusion  and  neat  conduction,  scattering 
and  random  walk,  and  wave  propagation.  In  this  paper  we  wish 
to  present  a  new  method  for  establishing  the  existence  of 
solutions  of  these  equations  in  the  cases  where  no  creation  of 
matter  is  involved.  This  method  is  based  upon  the  introduction 
of  a  new  class  of  physically  important  functions,  the  dissipation 
functions ,  Combining  these  new  functions  with  the  functions  we 
have  repeatedly  utilized  in  the  past,  the  reflection  and  trans¬ 
mission  functions,  we  easily  obtain  a  basic  conservation 
relation  upon  which  all  else  hinges.  The  uniqueness  proofs 
follow  conventional  lines. 
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DISSIPATION  FUNCTIONS  AND  INVARIANT  IMBEDDING — I 

Richard  Bellman  and  Robert  Kalaba 
Tne  RAND  Corporation 

0.  Milton  Wing 
Sandia  Corporation 

1.  Introduction.  In  a  aeries  of  papers  doting  from  1956  (see  1 
for  a  complete  bibliography),  we  have  used  the  theory  of  invari¬ 
ant  Imbedding  to  derive  a  variety  of  nonlinear  functional 
equations  for  the  description  of  radiative  transfer,  neutron 
transport,  diffusion  and  heat  conduction,  scattering  and  random 
walk,  and  wave  propagation.  In  thi**  paper  we  wish  to  present  a 
new  method  for  establishing  th-  ex ? sten  .  -  of  solutions  of  these 
equations  In  the  cases  where  no  creation  of  matter  is  Involved. 

This  method  Is  based  upon  the  introduction  of  a  new  class  of 
physically  important  functions,  the  dissipation  functions . 
Combining  ‘.h-rse  new  functions  with  the  functions  we  have  repeat¬ 
edly  utilized  In  the  past,  the  reflection  and  transmission 
functions,  we  easily  obtain  a  basic  conservation  relation  upon 
which  all  else  hinges.  The  uniqueness  proofs  follow  conventional 
i  1po° 

It  is,  of  course,  essential  for  the  rigorous  formulation  of 
the  Invariant  Imbedding  approach  to  present  these  existence 
thecrems.  The  results  we  present  are,  however,  of  deeper  signi¬ 
ficance  than  mere  rigor.  In  the  first  place,  upon  them  we  can 
four^  ~  approach  to  the  classical  linear  transput  piocesses, 

one  that  is  independent  of  characteristic  values  and  spectral 

2 

theory.  In  the  second  place,  as  Chandrasekhar  points  out, 
invariance  principles  cannot  always  be  guaranteed  to  yield  a 
unique  description  of  a  physical  process.  Consequently,  It  is 
essential  for  computational  purposes  to  know  that  we  have  a 
guaranteed  method  of  obtaining  the  physically  meaningful  solutions. 

2.  The  Mathematical  Model.  We  wish  to  consider  a  particle 
process  involving  N  types  of  particles  traversing  a  homogeneous 
"rod"  of  length  x.  We  suppose  that  when  a  particle  of  type  J 
penetrates  from  either  direction  a  segment  of  length  A,  then 
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incremental  particles  of  type 
opposite  direction. 


biJ  >0, 

the  process,  i.e.,  are  "absorbed,”  and 


fuA 


1  are  sent  back  in  the 

particles  disappear  from 

1,  A  additional 
ij 

particles  of  type  i,  i  /  j,  are  sent  through  in  the  same 

direction.  To  conserve  matter,  d .  <0  particles  of  type 

J 

J  are  added  to  the  original  stream,  where 
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JJ 


-  k 


'i  * 


‘u 


+ 1 


1.)  <  0. 

\J  \J  ) 


J  m  1,2, 


(1) 


Here  A  is  an  Infinitesimal  and  all  expressions  are  to  terms 
in  o  (A) . 

We  wish  to  determine  the  number  of  particles  in  state  i 
which  are  reflected,  transmitted  and  absorbed  per  unit  time, 
when  one  particle  of  type  J  per  unit  time  is  incident  upon  the 
rod.  Denote  these  quantities  by  r.  ,  (x),  t..(x),  and  >A.(x) 
respectively,  and  introduce  the  N  x  M  matrices  R{x)  -  (r^jC*)), 
T(X)  -  (tlj(x:;.  L(x)  -  (/^(x)),  B  -  (btJ ),  D  «  (d.p,  and 
the  diagonal  matrix  P  » 

Using  techniques  we  have  repeatedly  applied,  cf.  1 '  ,  we 
obtain  the  system  of  Riceati— type  differential  equations 


R 1 (x)  -  T ( x ) BT ( x ) ,  R(0 )  -  0, 

T ' ( x )  *  (D  +  R(x)B)T(x),  T(0 )  »  I,  (2) 

L'(x)  »  FT(x )  +  L(x)BT(x) ,  L(0)  »  0. 


It  is  this  system  of  equations  which  we  wish  to  discuss  for 
ail  x  >  0. 


Physical  Interpretation.  The  differentiation  of  particles 

can  be  made  on  the  basis  of  type,  energy,  wave  length,  angular 

direction,  and  so  on.  The  "rod”  can  represent  a  one-dimensional 
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rod,  a  slab,  or  generally  a  collection  of  strata,  as  in  .  ITie 
foregoing  is  thus  a  model  of  radiative  transfer,  multiple 
scattering,  and  similar  processes. 

4.  Conservation  and  Existence  Theorem.  Our  principal  result 


is  the  following. 
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Theorem.  Let 


blJ  ^  °’  dlJ  -  °’ 

i  /  j  > 

hi  >  ». 

(1) 

M(B  +  D  +  P)  -  0, 

(2) 

where  M  is  the  N  x  N  matrix  whose 

first  row  is 

composed  of 

ones  arid  whose  other  elements 

are  zeroes. 

Then : 

(a)  The  solution  of  (2.1)  exists  for  all  x  >  0  and  la 


unique . 

( b )  R ( x ) ,  T ( x )  and  L ( x )  are  nonnegative  for  all  x  >  0 . 

(c)  M ( R ( x )  4-  T(x)  +  L(x))  -  M  for  x  >  0.  {}) 

(d)  The  limits  R(a-),  ? (ec),  L(ou  )  exist . 

The  relation  in  ( ^ . c )  is  the  essential  conservation 
relation . 


jpj; _ Sketch  of  Proof.  It  la  easy  to  show  that  solutions  exist. 

in  a  small  x— interval,  0  <  x  <  xQ,  satisfying  the  conditions 
R(x),  T(x),  L(x)  2  °*  For  *?^aj.iple,  we  can  use  the  recurrence 
relations 


Rn+l(x)  "<^X  Tn(y)BTn(y)dy,  x  >  0, 

Tn+1  (*)  =  ^DX  eD^x“y^Rn(y)BTn(y)dy, 


(1) 


with 


R0(x) 


0, 


T0(x) 


I.  The  fact  that 


Tn+1<*>  ^  °  lf 

R  (y),  T n(y)  ?  0  follows  from  the  known  result  that  eu '  >  0 
if  d^  >0,  i  /  J ;  see  . 

If  m  now  establish  the  conservation  relation  of 
which  asserts  the  uniform  boundedness  of  r^j(x),  t^j(x), 
over  the  interval  [0,xj  ,  then  standard  techniques  in  the 
theory  of  nonlinear  differential  equations  will  penult  us  to 
assert  the  existence  of  the  solution  for  all  x  ^  0.  As 
mentioned  above,  uniqueness  fo.llf.ws  readily. 

We  have,  using  the  equations  of  (2.1), 

[M(R  +  T  +  L)]  '  -  M(R  ♦  T  ♦  L)BT  +  M(D  +  P)T.  (2) 


Consider  this  as  a  differential  equation  for  the  matrix 
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Z  «  M(R  +  T  +  L)  .  It  la  satisfied  over  [0,Xq] ,  the  interval 
of  existence  of  R(x),  T(x),  L(x),  by  the  constant  matrix 
Z  -  M  since  M(R (0)  +  T{0)  +  L(0))  -  MI  -  M  and 

0  *  M*  ■  MBT  +  M(D  +  P ;T  -  K(fi  +  D  +  F)T  -  0,  (5) 

upon  referring  to  (4.2),  we  see  that  (4.5c)  holds  over 
0  <  x  <  xQ. 

b.  Discussion.  In  this  paper,  we  have  applied  a  new  technique 
for  establishing  existence  of  solutions,  based  upon  conservation 
principles,  to  one  class  of  ffroblemn  .  In  forthcoming  papers,  we 
will  discuss  more  complex  problems  involving  time-dependence, 
particle— particle  interaction  and  variable  media.  We  shall  also 
discuss  uses  in  transmission  line  theory  and  wave  propagation. 
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